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effective single-particle equation? 
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The mapping, exact or approximate, of a many-body problem onto an effective single-body problem is one 
of the most widely used conceptual and computational tools of physics. Here, we propose and investigate the 
inverse map of effective approximate single-particle equations onto the corresponding many-particle system. 
This approach allows us to understand which interacting system a given single-particle approximation is actually 
describing, and how far this is from the original physical many-body system. We illustrate the resulting reverse 
engineering process by means of the Kohn-Sham equations of density-functional theory. In this application, our 
procedure sheds light on the non-locality of the density-potential mapping of DFT, and on the self-interaction 
error inherent in approximate density functionals. 
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One of the most widely used and successful approaches to 
many-particle physics is to map the many-body problem onto 
an effective single-body problem. Innumerable concepts and 
methods of theoretical physics derive from this general idea. 
The mean-field approximation, one of the most widely used 
approximation schemes in all fields of physics, is of this type, 
as well as the (in principle exact) mapping onto Kohn-Sham 
(KS) equations, used in density-functional theory (DFT). Here 
we introduce the concept of reverse engineering of single- 
particle equations [1], as a tool for discovering for which 
many-body system a given single-particle approximation be- 
comes exact. The idea is developed below in the framework 
of DFT, but the concept is completely general, and can be ex- 
tended to other single-particle methods. 

DFT [2] in the Kohn-Sham formulation [3] is a tool for cal- 
culating properties of many-body systems by means of a one- 
to-one mapping between the interacting system and a fictitious 
non-interacting one, with the same ground-state density n{r). 
All ground-state properties can in principle be expressed as 
functionals of this density. The success of DFT depends on 
the quality of the approximation to the exchange-correlation 
(xc) energy functional E,^c[n], which enters the KS equations 
through its functional derivative, the xc potential fxc(f )■ 
stems from the interactions in the original many-body system, 
and its functional form for Coulomb interacting systems, such 
as electrons in atoms, molecules, nanostructures and solids, 
is unknown. In order to construct viable approximations, a 
great deal of work has therefore been devoted to the deriva- 
tion of exact properties of the exchange-correlation functional 
and potential [4, 5] . The performance of a functional is judged 
by how close the density, and the observables calculated from 
it, are to that of the many-body system under study. 
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The idea of reverse engineering suggests a different mode 
of analysis, namely to ask: for what system does a given ap- 
proximation become exact? This question can be interpreted 
in two ways, one very common, the other being proposed here. 
To exemplify both, consider the local-density approximation 
(LDA). By construction, the LDA becomes exact for uniform 
densities. But in practice we rarely apply density function- 
als to uniform systems. Instead, we apply DFT to realistic 
inhomogeneous many-body systems, for which the LDA den- 
sity is r-dependent and approximate. The question we ask is: 
for which alternative many-body system is this approximate 
density the exact ground-state density? In this context, we 
call this alternative many-body system the interacting-LDA 
(i-LDA) system [6]. Reverse engineering of DFT refers to the 
process of constructing this alternative many-body system. 

The aim of this construction is not to simulate large sys- 
tems, for which nothing is gained by mapping a many-body 
system on another, equally complex, one. Rather, it shifts at- 
tention from the density predicted by the LDA to the external 
potential predicted by the LDA (i.e., the i-LDA potential), al- 
lowing us to investigate the LDA (or any other single-particle 
approximation to the many-particle Hamiltonian) in a novel 
and particularly detailed way. For example, we can now ask: 
how close is the corresponding i-system external potential to 
the true external potential? Which artificial features is the cho- 
sen approximation building into the system? As in other parts 
of science, reverse engineering enables one to understand the 
functionality and structure of the engineered device on a dif- 
ferent level, opening up new pathways for improvement. 

To illustrate the basic idea, we start from the LDA densities 
of the helium atom and of Hooke's atom, and invert the many- 
body Schrodinger equation to construct that external potential 
for which the LDA densities are exact ground-state densities. 
Comparison to the true external potentials of the helium atom 
(~ l/r) and Hooke's atom (^ r^) reveals the eiTors inher- 
ent in the approximate density and functional. By calculating 
explicitly the external potential for the i-LDA system, the pro- 
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posed method provides the exact i-LDA system Hamiltonian. 
This impHes that all properties of the i-LDA system (which 
themselves depend on the accuracy of the LDA) can, in prin- 
ciple, be directly calculated and compared to the ones of the 
exact many-body system. 

Studying relatively simple systems, such as heUum or 
Hooke's atoms, to understand approximations has proved 
fruitful in the past [7-9]. Earlier inversion schemes in DFT 
were used to find the Vx,c{r) that reproduces a given exact 
density, by inverting the single-particle Kohn-Sham equations 
[10]. This approach provides information on the exact ?;xc- 
The reverse engineering procedure, by contrast, aims at re- 
producing a given approximate density by inverting the many- 
body Schrodinger equation, providing information on the ap- 
proximate Wxc- Inversion of the many-body equation is a much 
harder task, which up to now had only been attempted for a 
one-dimensional model system within the adiabatic approxi- 
mation to time-dependent DFT [9]. To achieve this, we devel- 
oped inversion schemes that are substantially more accurate 
then previous ones (see discussion below). 

We consider the N-body Hamiltonian _ff = T -I- 14e + Vext, 
where T is the kinetic energy operator, Vee is the electron- 
electron interaction and t^xt = Vext{fi), is the external 
potential. By multiplying the Schrodinger equation H'ii = 
E^i from the left by and integrating over all but one of the 
coordinates we obtain 



**(ri...r^)(f + t>ee + t>ext) 



^{ri . . . rN)(rr2 ■ ■ ■ drvN = E 



N ' 



(1) 



From this we find, for the term concerning the external poten- 
tial. 



/ E «ext(ri)*dV2 . . . d'vN = «ext(ri)^ + 

{N-1)J **t;ex,(r2)*rfV2 . . . d^rN (2) 



n(vi) 2 

= Ven{ri)-j;^ + — I 7(ri,r2;ri,r2)wext(r2)'i^r2, 



(3) 



where7(ri,r2;ri,r2) = / ***(i^r3 . . . d^rAr. 

We combine these results to obtain an iterative relation for 
the external potential v^xtiri) that reproduces the target den- 
sity n'^s'^X^i). 
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-2 j 7j(ri,r2;ri,r2Kxt(»'2)c«V2]. 



To avoid having to calculate the integrals, we use the identity 

'^i(^l)^ext(»^i) - Eimiri) = 
-nJ^*{T + t>ee)*irfV2 . . . d^rN 



7i{ri,r2;ri,r2)vl^i{r2)d^r2, 



(5) 



and obtain the simple iterative relation 
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which is to be iterated with the Schrodinger equation (T + 
^ee + .V;xt)*^ = At convcrgencc, ni(ri) = n'^s'='(^i). 

and vl^i (ri) = Uext(^i) is the external potential that repro- 
duces this density via the many-body Schrodinger equation. 

We note that if ni(ri) is larger (smaller) than n'"S'^'(ri) 
then the potential must increase (decrease) at this point to 
bring ni-i-i(ri) closer to n'^''S'^'('ri). Therefore, iteration of 
Eq. (6) is expected to converge if Ei < 0. If Ei > 0, we 
replace Eq. (6) by 



(4) 



(^i) = ;r^^^["^(^i) - ""'"(^i)] + <t(r-i). (7) 

Scheme (6) converges relatively easily for the helium atom 
(where E <0) and scheme (7) for Hooke's atom (where E > 
0). We aid convergence in both cases by mixing v^^J^ with 
80% of Wg^,, and iterate until the relative error J d^r |nj(r) — 
n'"S'^'(r)|/ / d'^r n^^^^\r) has reached a desired level. 

At convergence, we obtain that external potential whose 
many-body ground state has the same density as was predicted 
by the approximate density functional for the true external po- 
tential. As this interacting system ground-state reproduces the 
LDA density then by the Hohenberg-Kohn theorem [11] this 
interacting system must be the exact interacting system cor- 
responding to the LDA. We check this property by indepen- 
dently solving the many-body Schrodinger equation with the 
converged alternative external potential, using a larger basis 
set than that used in the iterations. This additional consistency 
procedure allows our inversion scheme to be precise even in 
regions of space where the density is just lO"'^ a.u. (helium) 
and 10~^^ a.u. (Hooke's atom), i.e. orders of magnitudes 
smaller than previous inversion schemes, which, according to 
Ref. [9] attain an accuracy of up to 10~^ a.u. 

We now illustrate both inversion schemes, and the addi- 
tional consistency check, by applying our procedure to the he- 
lium atom and Hooke's atom, generating the approximate (tar- 
get) density from the LDA. For two electrons in a spherically 
symmetric potential the ground state can only be a function 
of the distance of each electron from the origin and the an- 
gle between the electron vectors. Hence we employ the basis 
(fiiji = Ri{ri)Rj{r2)V^rTTPi{cos{e))/{4iT). For Hooke's 
atom Ri (r) = Qi (r)e^"'' is a harmonic oscillator-like wave- 
function while for the helium atom it is a hydrogen-like wave- 
function Ri{r) = Qi{r)e~"^. The Qi{r) are polynomials of 
degree i created via the Gram-Schmidt procedure such that 
the Ri are orthonormal. 
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Hooke's atom is an interacting system of two electrons in 
the harmonic confining potential Uext = We use the 

LDA in the parametrization of Perdew and Wang [12] for v^c, 
and solve the KS equations for hio = 0.0365 Hartree. For 
comparison we also calculate the exact density from the exact 
interacting wave function using the method of Taut [13]. We 
then apply Eq. (7) to find the external potential of the inter- 
acting system that reproduces the LDA density, i.e. the i-LDA 
system [6]. 

In practice, an external potential found with a basis of 
e.g. 6^ functions may not reproduce the LDA density when 
we counter-check the scheme by solving the many-body 
Schrodinger equation of the i-LDA system using a larger basis 
of size e.g. 7^. Hence, we repeat the procedure with increas- 
ing basis size until we find an external potential that repro- 
duces the LDA density even with a larger basis. 

In Fig. 1 (a) the exact and LDA densities are plotted. The 
LDA and i-LDA densities are indistinguishable, but the differ- 
ence between the LDA and exact densities is large, and con- 
sequently the potentials of the true and the i-LDA system are 
rather different. Fig. 1(c) shows a substantial difference be- 
tween the i-LDA potential and the true external potential at 
large r, where the i-LDA potential grows more rapidly than 
the true potential. The LDA density is very different from the 
true density in this region (Fig. 1 (a) inset), with a relative er- 
ror (n(r) - nLDA(r))/„LDA _ gQ% r = 28 [14]. By con- 
trast, near the origin the i-LDA potential is slightly weaker 
than the true one. In between the two Umiting regions, the 
i-LDA potential and the true external potential cross various 
times. These crossings manifest themselves in a complex way 
in the behavior of the corresponding densities: Fig. 1 shows 
that a crossing in the potentials around r = 8 (Fig. 1 (c) inset) 
corresponds to a crossing in the densities (Fig. 1 (b)), whereas 
the crossings at r ~ 2 and r ~ 25 are not accompanied by a 
crossing in the densities. 

Next we consider the helium atom, for which we use the 
iterative scheme of Eq. (6) to calculate the i-LDA external po- 
tential. Here a basis of 7^ functions is required to find the 
external potential that satisfies our consistency check, i.e. re- 
produces the density when solving the Schrodinger equation 
with a larger basis (8"^ functions). 

As Fig. 2 shows, the LDA reproduces relatively well the 
exact radial probability density, but an underestimate at small 
r and overestimate at large r is evident. A closer view of 
large r (inset of Fig. 2(b)) shows that there is a larger discrep- 
ancy (49% on the range 4 < r/ao < 5.5 compared to 3.68% 
overall) between the LDA density and the exact density here. 
However we still find that the i-LDA system density is almost 
indistinguishable from the LDA density on this scale (0.78% 
error) and its overall error at 0.037% is comparable to that of 
0.0063% which we achieved for Hooke's atom. In Fig. 2(c) 
we again note that the external potential of the i-LDA system 
is substantially different from the true one for large values of 
r, where it is weaker than the true potential, causing the LDA 
density to be spread out shghtly more than the exact one. We 
also observe a crossing of both potentials and densities for 
r ~ 1.2 (Fig. 2(b)) where the radial probability density is 
high. 
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FIG. 1: Hooke's atom: huj = 0.0365 Hartree, basis size: S'^ func- 
tions (counter-checked with a 9'^ basis set, see text), (a) Comparison 
of LDA, exact and i-LDA densities. Inset: Zoom of the LDA and 
i-LDA densities at large r. (b) Comparison of LDA, exact and the 
i-LDA radial probability densities. Inset: Zoom of the tail of the 
density, (c) Comparison of the true external potential and the i-LDA 
external potential. Inset: Zoom of the potentials close to the origin. 



Both Hooke's and the helium atom thus exhibit a relation- 
ship between crossings in densities and crossings in poten- 
tials, but this relation is local only in regions where the ra- 
dial probability density is high. The systems considered here 
have spherical symmetry, hence the radial probability density 
47rn(r)r^ indicates the probabihty of finding a particle at dis- 
tance r from the origin. Thus, a single marked peak of the 
radial probability density at r indicates a high probability of 
finding a particle at that distance and a very low probability 
of finding the particle elsewhere. In this sense, the system is 
"almost classical" in the region of the peak. The charge den- 
sity itself does not have this property. The regions of high 
radial probability density can be interpreted as almost classi- 
cal, so that a local relationship between the 'particle position' 
and the external potential can be expected. 

By contrast, crossings in the potentials in regions with low 
values of the radial probabihty density do not necessarily cor- 
respond to crossings in the densities. These non-local regions 
are the ones in which the quantum nature of the system is 
more apparent. The identification of such regions is particu- 
larly interesting in view of recent investigations of semiclas- 
sical approximations to DFT [15]. The LDA density can even 
be higher than the exact density in regions where the i-LDA 
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FIG. 2: The helium atom: (a) Comparison of LDA, exact and i-LDA 
densities. Inset: zoom of the densities close to the origin, (b) Com- 
parison of LDA, exact and the i-LDA radial probabiUty densities. 
Inset: Zoom of the tail of the density, (c) Comparison of the true 
external potential and the i-LDA external potential. 

external potential is more repulsive than the true external po- 
tential. This dramatically highlights the non-local nature of 
the density-potential mapping in DFT 

Apart from the non-locality of the density-potential map- 
ping, the other main feature common to both systems is the 
substantial difference between the i-LDA and exact external 
potential at large r. For Hooke's atom in particular, the i- 
LDA potential diverges to infinity even faster than the exact 
potential [16]. We explain this as a consequence of the er- 
roneous asymptotic decay of the LDA exchange-correlation 
potential (exponential rather than 1/r), which itself is a con- 
sequence of the single-electron self-interaction error inherent 
in the LDA. Due to this error, the LDA xc potential is too weak 
in the asymptotic region, so that the i-LDA potential and ex- 
act potentials strongly differ in this region and, at least for the 
Hooke's atom [17], the i-LDA potential becomes much too 
confining at very large r's. However in this region n^'^'^(r) 
is still larger than the exact density. Hence the interplay be- 



tween many-body and single-particle self-interaction effects 
and between these effects and the non-locality of the mapping 
between density, Wxc and external potential is highly nontriv- 
ial. This is particularly evident for helium where the i-LDA 
potential is more attractive at the largest r depicted. Common 
to both systems is the drop in the i-LDA external potential af- 
ter the crossing at intermediate r's, this crossing is common to 
both potentials and densities. We speculate that this feature is 
caused by a less investigated consequence of self-interaction, 
causing the electron density to spread out more at smaller r in 
an attempt to minimize the self-Coulomb energy. As this starts 
to occur before the asymptotic region is reached, it involves 
many electrons, and can thus be interpreted as a consequence 
of the many-electron self-interaction error[18, 19]. 

The self-interaction problem in a single-particle framework 
has been studied for many years. Much more recently, it 
has become apparent that this single-particle self-interaction 
correction does not fully remove the self-interaction error in 
many-particle systems. The search for a clear indicator and a 
remedy for this problem has recently taken center stage as one 
of the main unsolved problems of DFT[18, 19]. Our reverse 
engineering prescription is a way to shed light on this com- 
plex problem, since the many-electron self-interaction error 
contributes to the difference between the true and the i-LDA 
potentials at intermediate distances. 

These considerations illustrate the concept and use of the i- 
LDA system. To construct this system we have developed two 
iterative schemes for inverting the many-body Schrodinger 
equation. Previous inversion schemes within DFT either were 
concerned with inverting the much simpler single-particle KS 
equation, or directed at a one-dimensional model many-body 
system. Our schemes attain orders of magnitude higher ac- 
curacy than previous schemes and are apphcable to three- 
dimensional many-body systems. 

Building on this technical advance, we have introduced the 
concept of reverse engineering in DFT and, more generally, 
in all many-body methods that introduce an effective single- 
body potential. This approach allows one to judge the per- 
formance and failures of, e.g., an approximate density func- 
tional or single-particle equation, by comparing two external 
potentials, thus revealing spatially resolved information on 
properties such as self-interaction errors and non-locality, in 
a physically transparent way. Another possible application of 
quantum reverse engineering is to design that external poten- 
tial that reproduces a desired density distribution in a given 
spatially inhomogeneous many-body system, which is an ex- 
citing prospect for the design of nano structured devices. 
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